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Abstract. We discuss the theory and experimental considerations of a quantum 
feedback scheme for producing deterministically reproducible spin squeezing. 
Continuous nondemolition atom number measurement from monitoring a 
probe field conditionally squeezes the sample. Simultaneous feedback of the 
measurement results controls the quantum state such that the squeezing becomes 
unconditional. We find that for very strong cavity coupling and a limited number 
of atoms, the theoretical squeezing approaches the Heisenberg limit. Strong 
squeezing will still be produced at weaker coupling and even in free space (thus 
presenting a simple experimental test for quantum feedback). The measurement 
and feedback can be stopped at any time, thereby freezing the sample with a 
desired amount of squeezing. 
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1. Introduction 

Squeezed spin systems |lj of atoms and ions have attracted considerable attention in 
recent years due to the potential for practical applications, such as in the fields of 
quantum information || and high precision spectroscopy ||. The basic principle is 
that the quantum correlations of squeezed spin states will outperform classical states 
in the same fashion as squeezed optical fields. Moreover, spin squeezing is related to 
the fundamental concept of entanglement and specifically represents many-particle 
entanglement |5|, || . There have been a number of proposals for spin squeezing and a 
variety of experimental results are now being observed. 

One area of research involves the use of broadband squeezed light , in which case 
the squeezed spin state is generated via quantum state transfer between nonclassical 
light and an atomic ensemble. This method has recently produced weakly squeezed 
states ||. In analogy with nonlinear optics, another proposal is the collisional 
interactions in a Bose-Einstein condensate (BEC). These represent a nonlinearity 
which will dynamically generate spin squeezing in the trapped state ||, [9) and also any 
out-coupled beams (lO). There are also schemes for direct coupling to the entangled 
state through intermediate states such as collective motional modes for ions jyj 
or molecular states for atoms [Q. A related proposal is the photodissociation of 
molecular condensates Jl3[ in analogy with the down-conversion process in quantum 
optics. There has also been experimental evidence that the ground state of a BEC 
confined in an optical lattice can be produced in an atom-number squeezed state p"4| . 

Production of spin squeezed states via quantum nondemolition (QND) detection 
has also been considered [[15[ and spin noise reduction using this method has been 
experimentally observed ]16| . QND measurements are also involved in the proposal 
Ji"7[ for the entanglement of two macroscopic atomic samples, which has recently been 
achieved |18| ]. These schemes represent conditional squeezing of the atomic ensembles. 
However, it would be preferable to develop unconditional, i.e., deterministically 
reproducible, squeezing especially in regard to the field of quantum information. The 
ideal would be to prepare the same squeezed spin state regardless of the measurement 
record. 

In this work, we develop the idea introduced in Ref. fl9| of achieving deterministic 
spin squeezing via quantum feedback, in analogy with optical squeezing via quantum 
feedback [^o| . Monitoring the output of a probe field sensitive to atoms in a particular 
internal state will conditionally squeeze the collective angular momentum of the 
sample. Unconditional squeezing is then obtained by using the measurement results to 
continuously drive the system into the desired, deterministic, squeezed spin state. We 
show (not done previously) that measurement schemes based in an optical cavity 
or free space both lead to the same system dynamics, i.e., that due to a QND 
measurement of one collective spin operator. 

It has been shown that a series of QND measurements followed by feedback 
can ensure perfect agreement between the number of atoms in two internal sates 
pi) . However this analysis does not take into account the quantum effects of the 
measurement back-action and feedback dynamics. The QND detection in our scheme 
will conditionally produce quantum correlations between the atoms, but it is not clear 
whether the feedback dynamics (introduced to produce unconditional squeezing) will 
adversely affect these correlations. In this paper, we perform a full quantum analysis 
to look at this question, and find the quantum limit to squeezing for the ideal situation 
of no atom losses. Further, we qualitatively explore how this limit is affected by losses 
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and what this implies experimentally. 

In Sec. 2 we briefly review the definition of squeezed spin states and the types 
of nonlinear interactions which produce them. The main body of this paper, Sec. 3, 
details the quantum dynamics of the continuous measurement and feedback scheme 
(3.1 and 3.2) and determines an experimentally feasible approach to implementing 
the required feedback. It also includes the effective nonlinear interactions produced 
by this scheme (3.3), details the numerical work (3.4), and a discussion of experimental 
considerations (3.5). Sec. 4 concludes and an appendix covers the results of inefficient 
measurement. 



2. Squeezed spin systems 

The collective properties of N two-level atoms are conveniently described by a spin- J 
system p2[ . This is a collection of 2 J = N spin-i particles where the internal states 
|1) and |2) of the kth atom represent the two states of the fcth spin-i particle. The 

collective angular momentum operators, J, are given by J a = J2k=i 3a {& = x i Di z )i 
where ja = 0a /2 and cr a are the Pauli operators for each particle. J obey the 
cyclic commutation relations [J x ,J y ] = ie X yzJz and the corresponding uncertainty 
relation (A J y ) 2 (A J z ) 2 > \\{J X )\ 2 . 

Equivalently, defining the operators a\ (i = 1, 2) for creating atoms in a particular 
internal state, we have J z = ^(a\a2 — a\a\), J+ = a\a\, J_ = a\a2 and the total 
number operator a\a\ + 0^02 equals a constant, i.e., N. J z represents half the 
total population difference and is a quantity which can be measured, for example 
by dispersive imaging techniques J x = (</+ + J-)/2 and J y = (J + — JJ)/2i 

represent the two quadrature-phase amplitudes of the collective dipole moment. 

For a coherent spin state (CSS) of a spin- J system, all the elementary spins are 
pointing in the same mean direction. Such a state satisfies the minimum uncertainty 
relationship with the variance of the two components normal to the mean direction 
equal to J/2. If quantum- mechanical correlations are introduced among the atoms it 
is possible to reduce the fluctuations in one direction at the expense of the other. This 
is the idea of a squeezed spin state (SSS) introduced by Kitagawa and Ueda (|], i.e. 
the spin system is squeezed when the variance of one spin component normal to the 
mean spin vector is smaller than the standard quantum limit (SQL) of J/2. 

There are many ways to characterize the degree of spin squeezing in a spin-J 
system. We will use the criteria of S0rensen and co-workers Q and Wang j2^] , where 
the squeezing parameter is given by 

2 iV(AJ ni ) 2 
4m (^ 2 > 2 + (J n3 ) 2 ' y 
where J n = n • J and n.j(i = 1, 2, 3) are orthogonal unit vectors. Systems with £ 2 < 1 
are spin squeezed in the direction n and it has also been shown that this indicates the 
atoms are in an entangled state [pip]. This parameter also has the appealing property 
that Q =£l = £ = l for a CSS J24J. 

It is well known that quadrature squeezing in optical or bosonic systems can be 
generated in a nonlinear Kerr (x^) medium. Kitagawa and Ueda Q identified similar 
nonlinear Hamiltonians that lead to squeezing in spin systems, namely 

Hi= X Jl (2) 

H 2 = X (Jl-J 2 -)/2i = x(J X Jy + JyJ X ), (3) 
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which can be seen as third-order nonlinear processes if we write J in terms of the 
atomic creation and annihilation operators. Hi can be thought of as twisting the 
z axis of the quasiprobability distribution. H 2 is then "two-axis countertwisting" , 
where the two axes bisecting the x and y axes are simultaneously twisted in opposite 
directions. Several recent studies have shown that Hi describes the nonlinear spin 
interactions in a two-component Bose condensate J5], |j| . 

As indicated in Ref. |Q, the variance of the squeezed spin component, and hence 
Eq. (Q), decreases to a minimum value corresponding to the maximally squeezed spin 
state before increasing again. The nonlinear interactions of Eqs. (Q) and (|j) twist 
the quasiprobability distribution of the spin state into, and then past, the optimally 
squeezed distribution. This is because the interactions can not simply be "turned off" 
when the optimal squeezing is reached. In the limit of large sample size, the minimum 
squeezing parameter for these interactions scale, respectively, as Hi : £ 2 cx N~ 2 ' 3 and 
H 2 :£ 2 (x N- 1 . 

In our proposal the required nonlinear interactions between atoms are introduced 
by a continuous measurement and feedback scheme. We not only find that the 
minimum squeezing parameter scales as N^ 1 (in the appropriate coupling regime), 
but that due to the nature of the scheme, the atomic sample could essentially be 
frozen in this maximally squeezed state (or less squeezed states) simply by ceasing the 
measurement and feedback. 

3. Quantum feedback scheme 

Let the two internal states of each atom, |1) and |2), be the degenerate magnetic 
sublevels of a J — | state, for example, the ground state of an alkali atom. For 
such atoms, transitions to the J = h excited states could then be employed for the 
QND measurement, as shown in Fig. tub). This figure shows an idealized experimental 
schematic for our QND measurement and feedback scheme (a), and the corresponding 
probe and driving transitions for each atom (6). The details of this scheme are 
presented in the following two sub-sections. 

We assume the atomic sample is prepared such that all the atoms are in one of 
their internal states and the temperature is very low (T w 0) such that we can separate 
out (and ignore) the spatial degrees of freedom. A fast 7r/2-pulse is then applied, 
coherently transferring all atoms into an equal superposition of the two internal states, 
which is an eigenstate of the J x operator with eigenvalue J. As described earlier, the 
CSS is a minimum uncertainty state so in this case variances of both J z and J y are 
J/2. 

3.1. Measurement 

The feedback scheme is based on a quantum nondcmolition measurement of the atom 
number difference between the states |1) and |2). The quantity to be measured is 
thus represented by the operator J z . Firstly, we consider placing the atomic sample 
in a strongly driven, heavily damped, optical cavity, as shown in Fig. 0(a). The 
cavity field is assumed to be far off resonance with respect to transitions probing state 
1 2), see Fig. [!](&). This dispersive interaction causes a phase shift of the cavity field 
proportional to the number of atoms in |2). Thus, a QND measurement of J z (since 
N is conserved) is effected by the homodyne detection of the light exiting the cavity 



Continuous QND feedback and unconditional spin squeezing 



5 




Figure 1. (a) Schematic experimental configuration. A cavity field (equally 
damped at rate n/2) of amplitude i\/3\ interacts with the atomic sample. The 
photocurrent I c (t) from the homodyne detection of the cavity output is combined 
with A(t) produced by a signal generator (SG). The combined signal controls the 
amplitude (AM) of an rf magnetic field that, together with a static field, drives 
J y . (b) Single atom diagram. A static B-ficld lifts the degeneracy of the magnetic 
sublevels. The far-detuned probe field monitors the collective population in state 
1 2) (and hence J z ). The radio frequency (rf) driving field, applied perpendicularly 
to the static field direction, induces magnetic dipole transitions between |1) and 
1 2) (thus driving J y ). 



In general, the interaction of a single atom with this far-detuned probe field is 
given by the Hamiltonian 

v = hn 2 ( Jz + ±)/4A (4) 

where fl is the Rabi frequency and A is the detuning. In the cavity configuration with 
iV atoms this becomes %x{Jz + N/2)b^b, where \ = g 2 /iA, and g is the one-photon 
Rabi frequency and b,tf are the cavity field operators normalized such that (b^b) 
is the mean photon number in the cavity. Note that we can choose an initial cavity 
detuning — Nx/2 to eliminate the N dependent term. For strong coherent driving we 
can use the semiclassical approximation b — > i\f3\ +6, where b now represents small 
quantum fluctuations around the classical amplitude The interaction is thus 

#cav =h X \P\Jz(-ib + itf), (5) 

where we have ignored the h X \f3\ 2 J z term by choosing an initial atomic detuning of 
-X\P\ 2 . 

We assume that the cavity field is damped through both mirrors at rate k/2, as 
shown in Fig. |l](a). The master equation for the combined atom-field system due to 
the interaction in Eq. (||) is thus given by 

Ptot = -x\P\[Jz(b - tf),p tot } + KT>[b]p tot , (6) 

where T>[r]p = rpr^ — (r^rp + pr^r)/2. Following the procedure of Sec. VII in 
Ref. pC)| we can adiabatically eliminate the cavity dynamics if the cavity decay rate 
K ^ x\P\AJ z , which requires k ^> x\P\VN (since initial AJ Z = a/ J/2). The evolution 
of the atomic system alone is then given by 

p = MV[J z ]p, (7) 

where M — 4x 2 |/3| 2 /k = Sx 2 P/n^K 2 and P = Huj\P\ 2 k/2 is the input local oscillator 
power and ui the frequency. The measurement strength M is equivalent to 2D in 
Eq. (22) of Ref. (25). Equation (0) represents decoherence of the atomic system due 
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to photon number fluctuations in the cavity field, with the result of increased noise in 
the spin components normal to J z . 

In the alternative case of a free space QND measurement, the analysis follows 
more closely that of Sec. IV A in Ref. ^6|. From Eq. (ji|), the interaction for many 
atoms in the free space limit is h9(J z + N/2)p*p. Here 9 is the phase shift of the 
probe field due to a single atom in state |2), and is given by 9 = ?iw7 2 /8AA/ sa t, where 
uj{— 2'kc/X) is the probe frequency, A the cross-sectional area, 7 the spontaneous 
emission rate and / S at = 2n 2 huj r y / A 2 for a two-level atom Note that the field 

annihilation operators p,p^ are normalised such that fiujp^p is the probe power, i.e., 
p'p is the photon flux operator with units of s _1 . The total interaction Hamiltonian 
(subtracting a phase shift of 9/2 for each atom) is in this case given by 

tffree = h9jjp. (8) 

The back-action on the atomic sample due to this interaction can be evaluated 
using the techniques of Sec. Ill B in Ref. (^8|. Assuming the probe field is a coherent 
state of amplitude i\g\, the system evolution is 

p=\ Q \ 2 V[e- i9J *]p~MV[J z ]p-i\Q\ 2 6[J z ,p\, (9) 

where the approximation requires 9AJ Z <C 1, which is equivalent to the requirement 
that 9\/N <C 1. The measurement strength is now M = \g\ 2 9 2 — P9 2 /hto, where 
P = fiui (jrp) is the mean power. The second term above simply causes a frequency 
shift in the atomic energy level difference, and can be ignored by choosing an initial 
atomic detuning of — \g\ 2 9. 

In both cases, the evolution of the system due to this type of interaction is given by 
the same expression MT>[J z ]p, i.e Eq. ([?]). The only difference is in the measurement 
strength, which is given by M = 8x 2 P/nu!K 2 for a cavity field and M = 9 2 P/Tiuj for 
free space. We can specifically look at the conditioned evolution of the system due to 
the homodyne detection of the probe field (this achieves the measurement of J z ). 

Continuing the analysis for the cavity regime as in p0[ , the conditioned density 
operator for the combined systems evolves as 

dptot = {-x\P\[J,Q>- b f ),Ptot] + KV[b]p tot }dt 

+ dW(t)^H[b]p tot , (10) 

where dW(t) is an infinitesimal Wiener increment defined by the Ito rules E[dW(t)] 
0: [cfW^)] 2 = dt, and H[r]p = rp + pr^ — Tr[(r + r^)p]p. We have also introduced 
a measurement efficiency ?y. For the setup in Fig. |lj(a), only one cavity mirror is 
monitored thus giving rj = 1/2 (assuming perfect photodetectors). In principle, 
however, one could monitor the output of both mirrors (using, e.g., a Faraday isolator 
at the input mirror) to obtain a perfect measurement, in which case n — 1. We will 
continue the analysis here with this assumption, and leave the treatment of n ^ 1 to 
the Appendix. 

Adiabatic elimination as before gives the stochastic master equation (SME) for 
the conditioned evolution of atomic system: 

dp c = MV[J z ]p c dt + \/MdW{t)H[J z ]p c . (11) 

In this equation p c is the density matrix for the atomic system conditioned on a 
particular realization of the homodyne signal given by 

I c (t)=2VM(J z ) c + C(t) (12) 
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where = dW(t)/dt is a white noise term satisfying E[£(t)£(t')] = 5(t — t') and 
E indicates the ensemble average. The average or unconditioned evolution is simply 
recovered by averaging over all possible realizations of the current, i.e. p(t) = E[p c (t)]. 
Note that these equations (|ll]) and (|l^) also apply for the perfect homodyne detection 
of the free space field, with the appropriate M. 

The dominant effects of the conditioned evolution ([ll]) on the initial CSS are a 
decrease in the uncertainty of J z (since we are measuring J z ) with corresponding noise 
increases in J y and J x , and a stochastic shift of the mean J z away from its initial value 
of zero. The latter shift can be calculated exactly (since d(J z ) c = Tr[J z dp c }): 



where the approximation assumes that (J z ) r ~ (such that I c {t) ~ dW/dt) which 
will apply in the next section. In effect, Eq. (O) is equivalent to a stochastic rotation 
of the mean spin about the y axis by an angle d<p ~ d ( J z ) / J. 

3.2. Feedback 

Monitoring the results of the measurement reduces the uncertainty in J z below the 
SQL of J/2 and increases the uncertainty in J y , but [as shown in Eq. (|l^)] the mean 
of J z stochastically varies from zero. The atomic system conditioned on a particular 
measurement result is thus a squeezed spin state, but the direction of the mean spin 
is stochastically determined. The unconditioned system evolution p = MV[J z ]p is 
obtained by averaging over all the possible conditioned states, and this leads to a spin 
state with (AJ Z ) 2 = J/2, i.e. the unmonitored measurement does not affect J z . In 
other words, the squeezed character of individual conditioned system states is lost in 
the ensemble average. This is illustrated by the states 2 and 3 in Fig. g. 

To retain the reduced fluctuations of J z in the average evolution we need a way of 
locking the conditioned mean spin direction. This can be achieved by feeding back the 
measurement results to continuously drive the system into the same squeezed state. 
The idea is to cancel the stochastic shift of ( J z ) c due to the measurement. This simply 
requires a rotation of the mean spin about the y axis equal and opposite to that caused 
by Eq. (|l3|). This is illustrated by state 4 in Fig. |. 

To make a rotation about the y axis proportional to the measured photocurrent 
I c (t), we require a Hamiltonian of the form 



where F(t) = X(t)J y /^/M and X(t) is the feedback strength. We have assumed 
instantaneous feedback because that is the form required to cancel Eq. (|l3|). Such 
a Hamiltonian can be effected, for example, by applying either the combination of 
a static and an amplitude-controlled rf magnetic field p9fl , as shown in Fig. |l|, or 
two-photon excitation by amplitude-controlled optical Raman fields. In either case, 
the amplitude of the driving field is controlled by the combined signal X(t)I c (t) [see 



d(J z ) c =2VMdW(t)(AJ z ) 2 c 
^2VMI c (t)dt(.J 2 z ) c , 



(13) 




(14) 




Continuous QND feedback and unconditional spin squeezing 




Figure 2. Schematic quasiprobability distributions for the spin state. The 
spin states are represented by ellipses on a sphere of radius J. The initial 
CSS, spin polarized in the x direction is given by state 1. State 2 is one 
particular conditioned spin state after a measurement of J z , while state 3 is the 
corresponding unconditioned state due to averaging over all possible conditioned 
states. The effect of the feedback is shown by state 4: A rotation about the y axis 
shifts the conditioned state 2 back to (Jz) c = 0. The ensemble average of these 
conditioned states will then be similar to state 4. 



where JCp = —i[F(t),p). The feedback terms of this conditioned evolution lead to a 
shift in the mean J z by an amount equal to 

d{J z )^ = -X(t)dW(t)(J x ) c /VM-X(t) 2 dt(J z ) c /2M 
- X(t)dt(J x J z + J z J x ) c 

m - X(t)I c (t)dt{J x ) c /VM. (16) 

As before the approximation assumes that (J z ) c = and also that the correlation 
between J z and J x is unchanged by the feedback. This is reasonable because it is 
initially zero due to the symmetry of the CSS and the conditioned states will remain 
symmetrical as shown in Fig. |[ 

Since the idea is to produce (J z ) c = via the feedback, the approximations above 
and in Eq. apply and we can find a feedback strength such that Eq. (|l^) is 
cancelled by Eq. @. The required feedback strength for our scheme is thus 

X(t) = 2M(j 2 z )J(J x ) c , (17) 

which is obviously dependent on conditioned averages. The feedback control, as 
expressed by Eqs. (|l4|) and ([17]), is essentially a form of state-estimation based feedback 
and has similarities to the method of Dohcrty et al |50|. Although it appears that 
the feedback depends directly on the instantaneous measurement current I c {t), the 
strength of this feedback is determined by the conditioned state expectation values. 
The current only appears directly in Hpo because of the assumptions we make about 
the conditioned state. 

Although Eq. ([l7|) is the ideal form of the feedback strength required to produce 
{Jz) c = 0, it may be not experimentally practical since it is dependent on conditioned 
expectation values. This requires classical processing of measurement signal in real 
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time to continually obtain the best estimates of the system variables J x and J 2 , which 
are then fedback to the system via Eq. (|l7|) and the Hamiltonian (|l4|) . What would be 
preferable is a series of data points or ideally an equation for \(t) that can be stored in 
a signal generator, as shown in Fig. |l](a), without having to do these side calculations 
during the feedback process. 

Applying the Markovian theory of Sec. IVA in Ref. we obtain the continuous 
QND measurement and feedback master equation for arbitrary A(i). It is given by 

f> = MV[J Z ] P - i\(t)[jy, J zP + P J Z ] + ^V[jy}p. (18) 

The terms in this equation describe, respectively, the measurement back-action, the 
feedback driving, and the noise introduced by the feedback. The master equation can 
also be rewritten in Lindblad form |E8|: 



—i 



c* Fit) + F(t)c, p] + V[c - iF{t)]p = Cp, (19) 



where the system operator c = \/MJ z and F{t) = X(t)J y / M as before. Note that 
these equations describe the exact unconditioned evolution of the atomic system where 
the feedback strength is arbitrarily defined by X(t). Equation ([lj]) thus describes one 
particular feedback scheme. 

To find an experimentally suitable expression for A(t) we first assume the feedback 
is successful in shifting the conditioned squeezed states to ( J z ) c = 0. This can be tested 
by looking at the purity, Tr[p 2 ], of the ensemble and conditioned states. Assuming 
perfect measurement, the conditioned spin state will remain in a pure state. Thus, 
if the ensemble average also has a high purity, it must comprise of nearly identical 
highly pure conditioned states. Solving the master equation, Eq. (|l8|), we find that 
the unconditional state does have a purity very close to one (see the numerical results 
section for details), and we are justified in applying the approximations (J x ) c — (Jx) 
and {J z ) c — (J z ) for the feedback strength. 

We can now proceed to find analytical expressions for (J x ) and { J 2 ) and hence 
X(t). To find an expression for (J x ) we note that its decrease from J is due to 
the increase of (AJ y ) 2 from J/2. The only contributions to this increase are due 
to the measurement, i.e., the first term in Eq. (|l8|), since the feedback operator 
(oc J y ) commutes with J y . Now, for (J) ~ O(J), the uncertainty in J y has not 
increased greatly and we can use the approximation (J x ) ss J (cos(A J y j J)) ss 
Jexp[— (J y ) /2</ 2 ]. Here the first approximation assumes that the angle defining 
the extent of the spin state in the xy plane is ~ AJ y /J, which requires J ^> 1, and 
the second approximation assumes near Gaussian statistics for J y . 

To find an expression for (J z ) we use assume that, for (J) ~ O(J), the atomic 
sample will approximately remain in a minimum uncertainty state. This is equivalent 
to assuming that the feedback, apart from maintaining ( J z ) = 0, does not significantly 
alter the decreased variance of J z that was produced by the measurement. This gives 
(J 2 ) ~ J 2 /4(j 2 }, where we have also used (J x ) ~ J. This last step is essentially 
a linear approximation represented by replacing J x with J in the commutator 
[J y ,J z ] = iJ x . From d t (j 2 ) = Tr[J 2 p] we obtain (J 2 ) w J 2 Mt + J/2, where the 
linear commutator has been used instead of the usual cyclic commutator. 

Substituting the approximation for (Jy) into the expressions for (J x ) and (J z ) 
we obtain 

(J 2 ) m (4Mt + 2/jy\ (20) 
{Jx) w Jexp(-Mt/2), (21) 
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and hence for the feedback strength we have 

\{t) « Me M */ 2 (l + NMt)-\ (22) 

So we now have a form of the feedback strength which can be applied experimentally. 
The question is - how well does it work? 

We can analytically approximate the degree of squeezing produced by the 
particular feedback scheme represented by Eq. fl22]) . For our model (ni,n 2 ,n 3 ) = 
(z,x,y) and the squeezing parameter, Eq. (|l]), becomes 



This leads to a minimum at t» sa 1/JW of 

£ in «e/AT, AT»1. (24) 

Thus, the minimum attainable squeezing parameter asymptotically approaches an 
inverse dependence on the sample size, i.e., the Heisenberg limit. This dependence is 
verified numerically in Sec. 3.4. 



3.3. Effective nonlinear interaction 

The quantum correlations of our scheme can be compared with those obtained via the 
different nonlinear Hamiltonians introduced ear lier, Eqs. (§ and (§). This requires 
the effective nonlinear Hamiltonian for our measurement and feedback scheme. The 
master equation in Lindblad form, Eq. (|l9|), is by definition separated into reversible 
(Hamiltonian) and irreversible (diffusive) terms. The nonlinear spin interactions 
produced in our scheme are thus given by the Hamiltonian 

H = \{c)F{t) + F(t)c) = \\{t){J z J v + J y J z ). (25) 

This is equivalent to the Hamiltonian H2 if we note that, for an initial CSS with 
mean spin in the z direction, Eq. (^|) produces a squeezed state in the xy plane Q. 
Thus, the appropriate "two- axis countertwisting" to squeeze the J z variance of of an 
initial CSS with mean spin in the x direction, has the form of Eq. (p5|). It is therefore 
not surprising that our scheme, since it has essentially the same nonlinear effect on 
the collective spin as H2, also gives a 1/N dependence for the minimum attainable 
squeezing parameter. 

The advantage of our scheme is that it does not rely on internal nonlinear 
spin interactions of the sample, such as those produced by collisions in a Bose- 
Einstein condensate. In addition to creating the squeezed state, these nonlinear 
interactions proceed to destroy the squeezing at later times. As opposed to internal 
spin interactions, which can not easily be "switched off" , the measurement laser and 
applied magnetic fields in our scheme can. This is simply represented by setting 
M = in our equations, which gives A(t) = and therefore the master equation 
becomes p — 0. Thus, no further evolution of the atomic system (in the rotating 
frame) occurs, and the spin is essentially "frozen" in the squeezed state. 

The exact timing for this procedure can be calculated prior to the run if 
all experimental parameters are known. For Heisenberg-limited squeezing the 
measurement and feedback should be applied continuously for a time t* = 1/M. On 
the other hand, if weaker squeezing is desired, the time is found by solving Eq. (p3|). 
For short times and N ^> 1 this becomes t ~ 1/MN£ 2 . In practice the limiting factor 
to squeezing will be the time for significant atom losses, as discussed in Sec. 3.5. In 
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Figure 3. Time dependence of the purity = Tr[p 2 ] (central curve), the 
expectation values (J x ) (A) and (j^) (B), and the squeezing parameter £ 2 (C) 
for N = 20. The solid lines represent the exact numerical solutions of the master 
equation with Eq. (E3) for the feedback strength. The dashed lines are the 
corresponding analytical expressions given by Eqs. @, Q and Q. 



any case, the degree of spin squeezing can simply be controlled by varying the duration 
of the nonlinear interaction. 

3.4- Numerical results 

In this section we will justify the approximations leading up to Eq. (|2^), as well as 
verifying the ./V" 1 dependence of Cmin- Regardless of the approximations, which were 
only used to determine the experimentally suitable form of A(i), the numerical results 
of this section are exact for the given feedback scheme. By this we mean that the 
master equation (|lj) [or Eq. ([l9])] is solved exactly for a feedback strength determined 
by Eq. ( pin . This can be done, for example, by using the Matlab quantum optics 
toolbox pl[ . 

The first approximation was that the conditioned averages could be replaced by 
ensemble averages. As stated earlier, this can be tested by looking at the purity of 
the unconditioned state due to the analytical X(t). By iteratively solving the master 
equation ( |T§| ) while updating X(t) at each time step [using Eq. (22)], we find the time 



evolution of the purity (= Tr[p 2 ]) and the expectation values ((A) = TrL4p]). In this 
way the exact solution for the squeezing parameter, 

e z =N{jl)/(J x )\ (26) 

is also found. 

The results of this simulation are shown in Fig. ^, where the purity is given by the 
central curve, the expectation values are given by curves A and B, and the squeezing 
parameter is given by curve C. We can clearly see that the purity remains near unity 
for the times of interest; at the numerical minimum of the squeezing parameter it 
equals 0.978. This implies that the measurement and feedback scheme has worked to 
produce nearly identical conditioned states. As discussed in Sec. 3.2, this justifies the 
replacement of the conditioned averages of Eq. ( |l7j ) with ensemble averages and also 
the further approximations to obtain Eq. (|22|). 

The analytical expressions for the expectation values, Eqs. ( pcf ) and (|2l]), and 
the squeezing parameter, Eq. (p3|), are also included in Fig. |3| (dashed lines). As 
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Figure 4. Typical time dependence of the squeezing parameter §^ for the 
continuous measurement and feedback scheme for three different sample sizes 
N = 20, N = 50, and N = 100. The dotted line represents no measurement and 
no feedback, i.e., the sample remains in the CSS. 
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Figure 5. Atom number dependence of the squeezing parameter minimum, 
-,-This figure presents the results of numerical solutions of the ME, and hence 
Eq. (ge|), with X(t) given by Eq. (0) (circles) and A(t) X 120% (diamonds). These 
values approach 3.33/N and 3.4971V respectively. The dotted line is the lower- 
bound for single-shot feedback also with a 20% error, while the analytical result 
for perfect continuous feedback e/N, Eq. (p4l), is the solid line. 



can be seen the analytical approximation (dashed curve A) for (J x ) is quite good, 
while the approximation (dashed curve B) for (J^) is only good for times up to the £| 
minimum. Since we are only interested in applying the scheme up until this point, this 
is acceptable. As a result though, the analytical minimum for the squeezing parameter 
is not a perfect fit to the exact numerical results. However, it is the correct order of 
magnitude and so we expect the analytical scaling predicted by Eq. (|2~i|), i.e. N -1 , to 
be correct. 

To reiterate, the analytical equation (|2^) for the feedback strength succeeds in 
maintaining the conditional squeezing (produced by the measurement) in a constant 
mean spin direction (hence the high values for the purity of the ensemble average 
state) . It is also of a suitable form for experimental realization. We now proceed to 
study the limit to the squeezing produced by this scheme. In Fig. || the squeezing 
parameter is plotted for three different sample sizes. The dotted line refers to the case 
of no measurement and no feedback. As can be seen, the squeezing becomes stronger 
as the number of atoms in the sample is increased. 
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To find how the minimum of £ 2 scales with increasing sample size, we extract £ 2 in 
numerically from the above simulation. The minimum attainable squeezing parameter 
approaches an inverse dependence on N as the sample size increases, as shown in Fig. ^|. 
In this figure, the analytical result for £ 2 lin , given by Eq. (24), is plotted along with the 
minima obtained from the numerical simulation. The analytical coefficient (e w 2.72) 
represents an error of ~ 18% compared to the numerical fit (3.33/A), which is to 
be expected from Fig. ||. Nevertheless, these errors only apply to scaling coefficients, 
not to the scalings themselves. Thus, for our continuous measurement and feedback 
scheme, we can conclude that the minimum attainable squeezing parameter scales as 
TV -1 and occurs after a time of the order of M -1 . The exact values for a given sample 
size can be obtained by numerical simulation. 



3.5. Experimental considerations 

Our continuous scheme, as opposed to a single-shot method, is very robust to any 
experimental errors in the feedback strength, as is shown in Figjq. The latter approach 
consists of a single (integrated) measurement pulse (see e.g., ]17[), followed by a single 
feedback pulse. If there is a relative error of e in the feedback strength, this will induce 
an error term (AJ| rr ) 2 ~ e 2 iV/4, which will dominate the total variance for N ^> 1. 
Thus £^ in will have a lower bound of e 2 , and will never scale as A -1 . On the other 
hand, as shown in Fig. | a large (20%) error in X(t) for continuous feedback does not 
affect the TV" 1 scaling. 

Although the analysis presented in this paper is based on instantaneous feedback 
(in order to determine analytically the limit to squeezing), including a finite delay 
time will have a limited effect as long as it is of the order of rfb ~ (NM) -1 . This is 
the fastest timescale in the equation (|22| ) for the feedback strength, and we discuss 
below the implications this time has for practical application. Furthermore, we have 
found the theoretical AT -1 scaling of the squeezing parameter, £ 2 , to be unaffected by 
inefficient measurements. This is detailed in the Appendix. 

Experimentally, the limit to squeezing will most likely be dominated by 
spontaneous atom losses due to absorption (and hence scattering) of QND probe 
light. For a detailed study of the effects of photon scattering in optical phase-shift 
measurements see the recent work by Bouchoule and M0lmer |3^] . Since the scattered 
photons carry information about the atoms that is not recorded, this process can also 
be treated in the same way as detector inefficiency (a discussion of this is presented 
in the Appendix, where we show similar scalings for £^ in to that of Ref. ]32|]). 

In the present discussion we take the more extreme approach of looking at atom 
losses from the spin-J system due to absorption of the probe light. In this way we 
estimate the limiting values of the spin squeezing parameter for both cavity and free 
space situations. From Eq. (|J) it can be seen that the single atom loss rate in the 
far-detuned limit will be T — 7f2 2 /4A 2 , where 7 is the spontaneous emission rate, Q is 
the Rabi frequency and A the detuning. For the cavity field interaction, fi 2 = g 2 |/3| 2 , 
M = &X 2 P / Tiu> K 2 and \ — g 2 /AA, and so T = KjM/g 2 . For the free space model, 
SI 2 = 7 2 //2/ sat , P = I A, and M = P9 2 /huj where 6 and 7 sa t are defined before 
Eq. (||), and so we obtain T — 16ir 2 AM/\ 2 . The single atom loss rate can therefore 
be expressed in general as 

T = aM, (27) 

where a = K^/g 2 for the cavity and 16tt 2 A/X 2 for free space, and so the number of 
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atoms lost by time t can be estimated by 

AN = TNt = aMNt. (28) 

Therefore, by the time (~ 1/M) to reach the Heisenberg limit we have lost AiV ~ aN 
atoms. From Sec. 3.3, the time to reach an arbitrary amount of squeezing is 
t ~ 1/NM£ 2 , so we can say that the total number of atoms lost scales as AN ~ a/£ 2 . 

Determining the maximum number of atoms that can be lost before destroying 
a given value of squeezing will thus give us some realistic constraints that should be 
placed on a. Recently, Andre and Lukin analyzed the squeezing produced by the 
counter-twisting Hamiltonian, Eq. (^) , including the effects of dissipation |34) . They 
showed that for a given amount of spin squeezing £ 2 = 1/s, where 1 < s < N, the 
maximum number of atoms that can be lost without destroying the squeezing scales 
as AN max ~ (N/s) Ins Q. In other words, to achieve the desired level of squeezing, 
we need to set AN ~ A7V max . This becomes a/£ 2 ~ A^ 2 ln(l/£ 2 ). Given an arbitrary 
a, the amount of squeezing produced including atom losses therefore scales as 

~ ^c7/N. (29) 

So to produce any squeezing we need a < N, for £ 2 ~ 1/y/N we need a ~ 1, and for 
the Heisenberg limit we need a ~ 1/N. 

For the case of the cavity field interaction, a — / g 2 , so the requirement for 
Heisenberg limited squeezing becomes g 2 ~ Nkj. This can only be fulfilled in the very 
strong coupling regime, which, although difficult, has been achieved experimentally 
[ ^3| . Looking at the other regimes, we find that £ 2 ~ 1/y/N squeezing will be observed 
for g 2 ~ K7, and the generation of any squeezing, i.e., £ < 1, will be observed for 
g 2 N > K'y. These are the same results as found in Ref. |54|. As might be expected, 
the free space model will not reach the Heisenberg limit. In this case a = 16ir 2 A/A 2 , 
so the requirement becomes A ~ A 2 /16ir 2 N which is not possible. The limiting 
value for squeezing in this case would be near 1 / v^V, since this would be produced if 
A ~ A 2 /167r 2 , and to achieve any squeezing we only require A < NX 2 /16tt 2 . So for 
the free space regime we still expect a significant amount of squeezing to be possible. 

We can also determine the constraints that need to be placed on the probe laser 
by recognizing that our analysis is based on a far-detuned optical probe with A 3> J7. 
This implies that T <C 7, and thus from Eq. (|27|), M <C "f/a. The measurement 
strength for both schemes can be re-expressed as 

M = j 2 P/huA 2 a 2 , (30) 

and hence our laser requirements are 

P/A 2 < huja/j. (31) 

Take for example N = 10 7 Cesium atoms, where A = 852 nm (uj/2ir = 352 THz) 
and 7 = 5 MHz. Then to obtain spin squeezing in either the free space or cavity 
regimes, we only need P <C 10 _18 A 2 . On the other hand, to reach Heisenberg limited 
squeezing (a ~ 1/^V) we have the added requirement that P -C 10~ 33 A 2 as well as 
strong cavity coupling. Looking at a detuning of A = 1 GHz, this restricts the laser 
power to femtoWatts, which should be possible [p5[ . 

Finally, realistic feedback delay times will give some additional restrictions on the 
possible squeezing. As mentioned earlier, the feedback needs to act on a timescale of 
order 1/NM, which becomes tr, ~ TiujA 2 a 2 /NPj 2 using the general expression for M. 
For P = 1 fW and A = 1 GHz, this time is of the order 10a 2 /A s. Typical feedback 
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delay times are of order 10~ 7 s j35j, so unless the number of atoms is restricted to a 
few hundred, the Heisenberg limit (where a ~ 1/N) will not be feasible. 

During the preparation of this paper we became aware of an analysis of the 
entanglement in bad cavities by S0rensen and Maimer (3(|. Although the scheme 
described in that work is very different from ours (in fact it is similar to Ref. it 
also reaches similar results. All three schemes require very strong cavity couplings to 
produce Heisenberg limited squeezing, but weaker (and perhaps more useful) squeezing 
can be produced in bad cavities and free space. 

4. Conclusion 

In conclusion we have detailed a scheme for achieving, and controlling, the spin 
squeezing of an atomic sample by a continuous quantum feedback mechanism using 
an optical probe. This scheme has the advantage over previous schemes involving 
QND measurement |l^, |l6| in that it produces unconditional, or deterministically 
reproducible, squeezing. This means that the same squeezed spin state, which is also 
an entangled state (5[ will be produced regardless of the particular measurement 
record. It is therefore a more practical state for use in quantum information. 

The theoretical squeezing (or entanglement) parameter approaches the Heisenberg 
limit of £ 2 ~ 1/N for very strong cavity couplings, a limited number of atoms, 
very large detuning and small laser power. This 1/N scaling indicates a stronger 
squeezing mechanism than collisional interactions in a Bose-Einstein condensate where 
the scaling is TV" 2 / 3 §, g|. By ceasing the measurement and feedback interactions 
when the desired squeezing parameter has been reached, the squeezed state could be 
maintained indefinitely. Thus, once this state is generated, there will be no further 
decoherence due to the interactions of its production. 

Significant squeezing will still be produced at weaker cavity couplings, and even 
in free space, with correspondingly less stringent experimental requirements. The free 
space QND measurement regime therefore presents a relatively simple, and interesting, 
experimental implementation of quantum feedback, an area which is still in its infancy. 
Also, the weaker the squeezing the less vulnerable the state is to further decoherence 
after its production. The Heisenberg-limited squeezed state will be severely affected 
by the loss of a few atoms, whereas less strongly squeezed states can lose many without 
serious damage [Q. The states produced in a free space scheme might therefore be 
adequate for practical quantum information, where other experimental factors could 
destroy the more sensitive states. Then, if necessary, one can devise entanglement 
purification procedures for atomic samples |S7| . Finally, achieving and controlling a 
spin squeezed state could also be useful in obtaining tomographic phase space pictures 
of nonclassical states of atomic samples . 

Appendix A. Inefficient measurements 

In real experiments there is a finite efficiency involved in the measurement process. 
We now include this in our analysis to see whether the degree of spin squeezing is 
significantly affected. Basically, only a certain proportion (n) of the probe beam is 
actually detected and then used for feedback, but the undetected proportion (1 — 77) 
still causes back-action on the system. This could be due to only monitoring the 
output of one side of a cavity [as in Fig. 0(a)], or imperfect photodetectors, or even 
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undetected photons scattered off the atomic sample (which is treated more rigorously 
in Ref. §§). 

Let us first look at the effect of an inefficient measurement of J z , before adding the 
feedback. The stochastic master equation dp c is now made up of two contributions: 
the actual homodyne detection of a proportion n of the probe light 

rfD[c} Pc dt + y/rjdW(t)H[c]p Ci (A.l) 
plus the back-action due to the probe light that wasn't detected 

(1 - n)V[c]p c dt, (A.2) 
where c = \[M.J Z as before. Similarly, the homodyne photocurrent is now given by 

I c (t)=r,(c + J) c + ^im (A.3) 

The stochastic shift in (J 2 ) c due to an inefficient measurement is therefore 
2y/i]MdW(t)(A J z )c- If we again assume (J z ) c ~ (which will be true when the 
feedback is included) then I c (t) ss ^/rjdW/dt and so we again obtain the approximation 
in Eq. £L3|). 

Including feedback to cancel this stochastic shift, we obtain the same feedback 
Hamiltonian ( |l4| ) as before, although note that the current I c (t) is now defined by 
Eq. ( |A.3| ). Thus, the shift in the mean J z due to feedback is still given by Eq. ( fl6| ) 
and the optimal feedback strength is therefore also unchanged from the form given in 
Eq. ([IT]). It is the assumptions made for obtaining the analytical expression for A(i) 
that are now complicated. 

The master equation for the inefficient measurement and feedback is given by 

p=(l-T) + r,)V[c]p - i[F, cp + pJ] + -V[F]p, (A.4) 

V 



where c = \/MJ z and F — \(t)J y /\J M as before. Thus, the total back-action due to 
the probe does not depend on the detection efficiency since it is made up of both the 
detected (77) and undetected (1 — 77) contributions. The Lindblad form of the master 
equation, Eq. (|l9|), now also has the additional term [(1 — rj)/rj\T>[F\p. 

To simplify the form of X(t) in Eq. (|l^) we again assume the feedback is 
successful and replace the conditioned averages with ensemble averages. To find ( J2) 
we earlier assumed a near minimum uncertainty state, but this cannot be true for 
imperfect measurements. To use (J z ) ~ J 2 /4(j^) for 77 < 1 we must therefore 
calculate the average Jy only using the evolution due to the actual measurement and 
feedback, i.e. we ignore (1 — rj)T>[c\p. This gives (Jy) sw J 2 r]Mt + J/2 and then 
(J z ) w {ArjMt + 2/J)~ 1 . To find (J x ) we earlier noted that it decreased indirectly 
due to the measurement back-action. As shown above, the total back-action due to 
the probe beam is independent of detection efficiency and so we can keep the earlier 
approximation ( |2l| ) for (J x ). 

These approximations for the averages this time lead to a feedback strength of 

A(i) w Me M */ 2 (l + rjNMt)- 1 . (A.5) 
The approximate squeezing parameter then becomes 

^^e Mt (l + riNMt)- 1 , (A.6) 
which has a minimum at t» w 1 /r/M of 

e^^e^/N, AT»1. (A.7) 
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So for rj < 1, the minimum attainable squeezing parameter is greater than that for 
rj = 1. However, finite efficiency only affects the coefficient of this minimum - the 
iV _1 dependence is unaffected. Note that this is the same scaling with the number 
of atoms as found in Eq. (59) of Ref . |52| , which confirms that photon scattering can 
also be regarded as a measurement inefficiency. 
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